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Abstract 

We consider the modified Korteveg de Vriez equation on the whole line. Initial data is 
real and step-like, i.e. q(x, 0) = for x > and q(x, 0) = c for x < 0, where c is arbitrary 
real number. The goal of this paper is to study the asymptotic behavior of the initial- 
value problem's solution by means of the asymptotic behavior of the some Riemann-Hilbert 
problem. In this paper we show that the solution of this problem has different asymptotic 
behavior in different regions. In the region x < —6c 2 t and x > 4c 2 t the solution is tend to 
c and correspondingly. In the region —6c 2 t < x < 4c 2 t the solution takes the form of a 
modulated elliptic wave. 

1 Introduction 

Initial value problems with step-like initial function have very long story beginning from the papers 
by A.V. Gurevich, L.P. Pitaevsky [5] and E.Ya.Khruslov [S] in a middle of 70th. There are many 
papers devoted to different aspects of these problems. For the present time there are a few full 

and rigorous results on an asymptotic behavior of such problems [10], [?]-[?], [IS], ,. We pay 

our attention here to a problem which was not considered elsewhere and give results in a rigorous 
form using the method of the matrix Riemann-Hilbert problem and the steepest descent method 
for oscillatory matrix RH problems. We consider the modified Korteweg de Vries equation on 
the whole line. Initial datum is a step-like, i.e. q(x, 0) = for x > and q(x, 0) = c for x < 0, 
where c is an arbitrary real number. Without loss of generality we put c > 0. This problem 
can be considered as a shock problem. The goal of this paper is to study asymptotic behavior 
of the Riemann -Hilbert problem whose solution gives the solution of the initial-value problem. 
In this paper we show that the solution of the shock problem has different asymptotic behavior 
in different regions. In the regions — oo < x < —6c 2 t and 4c 2 t < x < oo, the solution is trivial, 
i.e. it is equal to c and respectively. In the region — 6c 2 t < x < —4c 2 t the solution takes the 
form of a modulated elliptic wave of finite amplitude. Thus for a large time the solution has finite 
amplitude in the first two regions while in the third region (4c 2 t < x < oo) it takes the form 
of a vanishing (as 0(t -1 / 2 )) self-similar wave. The development of the Riemann-Hilbert method 
for the shock problems arising for integrable PDEs on the whole line with the different finite-gap 
boundary conditions as x — > ±oo goes back to the works done in 80-90s by R. Bikbaev, P. Deift, 
V. Novokshenov, and S. Venakides. All those results was devoted to the initial- value problem 
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with self-adjoint Lax operators. Most recently, an implementation of the RH scheme to the shock 
problem and the evaluation of the long-time asymptotics of the solution to the focusing nonlinear 
Schrodinger equation on the whole line, where the Lax operator is not self-adjoint and the initial 
function was chosen in the form q(x,0) = Aexp(i[i\x\), was done in [16J. It is worth mentioning 
that our shock problem is different from that in [16] as well as our construction of the phase g - 
function. 

The inverse scattering transform method (1ST) for solving initial-value problems for nonlinear 
evolutionary equations, discovered in 1967 [TJ, turned out to be a very powerful tool, which 
allowed to obtain a huge number of very interesting results in different areas of mathematics and 
physics. At the beginnig of 90th a new great achievement in the further development of the 1ST 
method have been done by P.Deift and X.Zhou. It is a nonlinear steepest descent method for 
oscillatory matrix Riemann-Hilbert problem. With the new method it came a nice possibility 
to rewrite known asymptotic results for different nonlinear integrable models in the rigorous and 
transparent form (sf . [TTJ , [12] , [32] ) and obtain numerous new significant results in the theory of 
completely integrable nonlinear equations, random matrix models and orthogonal polynomials, 
integrable statistical mechanics. Our goal is to bring new results in theory of the shock problems, 
especially in the case of non self-adjoint Lax operators, and some development of ideas, given 
recently in [T7], [TS] in the direction of the strengthening of the nonlinear steepest descent method 
for oscillatory matrix Riemann-Hilbert problem. 

Let us consider the problem 

q t + 6q 2 q + q xxx = (1.1) 
I 0, x — > +oo 

q(x,0) = q (x) ^ { (1.2) 
I c, x — oo, 

where qo(x) is arbitrary step- like function, including the discontinuous case: qo(x) = c for x < 
and qo(x) = for x > 0. We suppose that the solution q(x, t) of this problem exists for 
K+. To study the initial value problem (ll.ip - (11.2p we will use the Lax representation of the MKdV 
equation in the form of over-determined system of differential equations: 

$ x + ik<r 3 $ = Q(x,t)$ (1.3) 

$t + Uk 3 <T 3 Q = Q(x,t,k)Q, (1.4) 
where $ = <&(x,t, k) is a 2 x 2 matrix-valued function, 

°> ■= (J -°i) ■ w-=U*, t) Vy < L8 > 

Q(x, t, k) = 4k 2 Q(x, t, k) - 2ik(Q 2 (x, t, k) + Q x (x, t, k))a 3 + 2Q 3 (x, t, k) - Q xx (x, t, k), (1.6) 

and k G C. The equations ( 11. 3p and ( ll.4j) are compatible if and only if the function q(x,t) satisfy 
the MKdV equation (II. ip . To apply the inverse scattering transform to the problem ( II. ip . (II. 2p 
we have to construct matrix valued solution of this equations defined by their asymptotics: 

$(x,a) = e- (M3l,ff3 + 0^j, x^+oo, lmfc = (1.7) 
V(x,t,k) = E(x,t,k) + (- J , x^-oo, ImA; = 0. (1.8) 
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Here E(x,t, k) is the solution of the linear differential equations 

E x + ikcr 3 E = Q C E (1.9) 
E t + 4ik 3 (T 3 E = Q c (k)E, (1.10) 

where 

Qc ■= ( °„ nl (1-H) 



c 0, 

Q c {k) =4k 2 Q c -2tkQ 2 c a 3 + 2Q 3 c . (1.12) 
We chose the solution E(x,t, k) as follows: 



E(x,t,k) = - 



( m + *(*) - \ 

x(re) x(k) 

V x(A;) -xM x(A;) + xM/ 



where 



k — zc 



X(fc) = v^ 2 + c 2 , Q(k) = 2{2k 2 - c 2 )X(k), x{k) = \\ - — - (1.14) 

X{k) and x{k) are analytic in the complex plane cut along the segment [ic,—ic], i.e. k G C \ 
[— ic, ic], and branches of roots are as follows: X(l) > 0, x(oo) = 1. 
The solution $(x, t, k) can be represented in the form: 

oo 

®(x,t,k) = (e~ ikxa3 + J K( y x,y,t)e- ikyas dyy- 4ik3t ' 73 , (1.15) 

X 

where the kernel K(x,y,t) is chosen to be so that the first factor satisfies the x— equation (11.31) 
for all t, and the second factor satisfies the t— equation (11.41) for x = oo. Then, by the same way 
as in [15J, we prove that $(x,t, k) satisfies both equations (11.31) and (11.41) . The existence of the 
solution represented by the transformation operators with the kernel K(x,y,t) is proved in |15j . 
By the same manner another solution takes the form: 

X 

V{x,t,k) = E(x,t,k) + J L(x,y,t)E(y,t,k)dy (1.16) 

— oo 

with some matrix kernel L(x, y, t). Omitting rutin details of the proof of these representations we 
formulate below properties of the solutions. 

The matrices Q(x,t, k) and ^(x,t, k) defined by (11. 7p and (II. 8p and their columns $j(x, t, k) 
and ^j(x, t, k), j — 1, 2 have the following properties: 

1. determinants: 

det $(x, t, k) = 1, det ^(x, t,k) = l. 

2. analyticity: 

$i(x, t, k) is analytic in k G C_, $2(x,t, k) is analytic in k G C+, 

^i(x,t, k) is analytic in k G C + \ [0, ic], ^/ 2 (x,t, k) is analytic in k G C_ \ [— ic, 0], 

^/i(x,t, .) and \J/ 2 (^c, ^ , •) have continuous extensions to (— zc, zc)_ U (— ic, ic) + . 
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3. symmetries: 



gggQMj) =*n(x,t,k), ^ 22 {x,t,-k) =V n (x,t,k), 
^i 2 (x,t,fc) =-* 21 (x,t,A;), *i 2 (x,t,-A;) = -* 2 i(x,t, k), 

$i2(x,t,k) =-$2i(x,t,k), ®i 2 (x,t, —k) =-§2i(x,t,k), 
t, -fc) = <f>ji(x, t, k), j, I = 172. 

4. large fc asymptotics: 

}-i + og). — , r m „o ; 

Xlr^T- f l.\ p +ikx+iik z t 1 /1\ 

«•< r i+o U> lmi - 0; 

5. jump: 

»_(*,*,*) = *+(.,U) ( J J ) , k e (-ic,ic), where » ± (,, t ,t) are the boundary val- 

ues of the matrix *&(x, t, k) from the left (+) and from the right (— ) of the oriented downwards 
interval (—ic, ic). 

The matrices t, k) and *&(x, t, k) are the solution of equations (ll.3p and ( 11. 4p . Hence they 
are linear dependent, i.e. there exists the independent on x,t matrix: 

T(k) = k)^(x,t, k), fcGl (1.17) 

which is defined for real k. Some of elements of this matrix have a larger domain of the definition. 
Indeed, using (11.17P we find 

T\i{k) = det( 1 i r i, $ 2 ) (1.18) 

T 21 (fc) = det($ 1 ,tf 1 ) (1.19) 

Ti 2 (fc) =det(^ 2 ,$ 2 ) (1.20) 

T 22 (A;) = det($i,* 2 ). (1.21) 
Then the above properties of the solutions $>(x,t, k) and ty(x,t, k) yield: 

• Tn(fc) is analytic in k G C+\[0,ic] and has a continuous extension to (0, ic)- \J(0,ic) + , 

• T 22 (k) is analytic in G C_\[0,zc] and has a continuous extension to (— zc, 0)_ |J(— zc, 0) + , 

• T 2 i(/c) is continuous in G (—00, 0) |J(0, — ic]- |J[— ic, 0) + lj(0, +00) 

• Tu(k) is continuous in k G (— 00, 0) |J(0, ic]- {J[ic, 0) + IJ(0> +°°) 
and 



T 2 2(k) = T 11 (k), T 22 (-k)=T n (k), 
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T 12 (k) = -T 21 (fc), T 12 {-k) = -T 21 (k), 



a(k) = T u (k), 
b(k)=T 21 (k). 



• T jk {-k) = T jk {k), j, k = 1,2 
Denote 

Define the reflection coefficient 
It has the following property: 

r(—k) = r(k) 

The columns of the matrices $ and ^ have the following extra properties: 
(*i)_(x,t,A;) (*i)+(x,t,A;) 



r(/c) 



a(/c) 



6. 



7. 



a-(k) a +(k) 
(* 2 )-(x,t,k) (V 2 )+(x,t,k) 



a_(fc) 



a + (/c) 



where 



/i(/O$ 2 0M,A0> fce(0,ic) 
f 2 {k)$ 2 {x,t,k), fce(0,-ic) 

fe g (0,zc) 



a-(k)a + (k) 



fi(k) 



a_(fc)a + (/c) 



-fi(k), ke(-ic,0) 



2 The Basic Riemann— Hilbert problem 



The scattering relations (11.17P between matrix-valued functions fc) and fe), and also 

extra properties 6, 7 can be rewritten in terms of the Riemann-Hilbert problem. To do so, let us 
define matrix- valued function M (£, t, k) by putting 



( f *i(x,t, k) 
a(k) 



e <tf(fc,O j $ 2 ( X; t? fc) e -^(*>0 J , jfc G C+\[0, ic] 



Af(£,t,A;) = 4 



a(k) 



(2.1) 



fc G C_\[-ic,0], 



where x = 12ft and 0(fc,f) = 4fc 3 + 12&f (f = x/12t). To make the paper more transparent we 
consider below only the shock problem when the initial datum is discontinuous: 



q (x) 



0, x > 
c, x < 0. 
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Then 



r(k) 



x 2 {k) - 1 
! (fc) + l' 



(2.2) 



where is defined by (11.141) . are analytic in G C\[— zc, zc]. The transition coefficient a^ 1 {k) 
is bounded in G C + \[zc, 0] because the function a(k) is equal to zero nowhere. In this case we 
have: 

/ 2 (A:) = f 1 (k)=r_(k)-r + (k), k G (-ic,ic) (2.3) 



IC 







-ic 



Re k 



Figure 1: The oriented contour E 

Let us define the oriented contour E = R U [ic, —ic] as on the figure 1. Then the matrix 
M(£, t, k) solves the next Riemann-Hilbert problem: 

• matrix valued function M(£, t, k) is analytic in the domain C \ E; 

• M(£, t, k) is bounded in neighborhoods of the end points ic and — ic; 

• M_(£, t, k) = Af + (f , t, k)J(£, t,k), k G E = R U [ic, -ic], where 



1 _ r (£;)e- 2i * e ^) \ 

_ r ( k yit9(k,z) i + | r (jfe)|2 j 



k G R\{0} 



J(£,t,k) = { 



\ 



( 1 f(k)e- 2ite ^ \ 



k G [0,zc] 



I \ 1 

M(£,t,k) = I + Oik- 1 ), k^oo; 



k G [0, 



(2.4) 
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where r(k) = b(k)/a(k), k G R is given in (12. 2p . and 



f(k) := fi(k) = f 2 (k) = — -— — — - = = = , k G [ic, -ic\. 

a_{k)a + {k) Q _( A .) a+ (Q 

If the initial datum is a generic step-like function then a(k) can have zeroes in the domain of 
analyticity. In this case the matrix M(£, t, k) will be meromorphic and residue relations between 
columns of the matrix M(£, i, k) must be added. 

Now we forget about the origin of the Riemann-Hilbert problem and suppose that the oriented 
contour £ and functions ( 12. 2 p are given. The following theorem take place. 

Theorem 2.1. Let the oriented contour £ and functions \2. fy) with 



^ = v hr- 

V k + ic 

be given. Then the Riemann-Hilbert problem has a unique solution M(£,t,k). The function q(x,t), 
given by the equations, 

q(x, t) = 2i lim [kM(x/12t, t, fc)] 12 , (2.5) 

fc— >oo 

satisfy the MKdV equation ( I J. 1\) and the initial condition 



q(x,0) = q (x) 



x>0 
x < 0. 



The proof of this theorem almost the same as in [T3], if we takes into account that given 
functions a(k), b(k) via the function x(k) are in the one-to-one correspondence with the shock 
function qo(x). 



3 Long time asymptotic analysis of the Riemann-Hilbert 
problem 

3.1 Steadiness region x < —6c 2 t 

The jump matrices J(£,t, k) depend on exp{±2it#(/c, £)}. Hence the signature table of the imag- 
inary part of 6(k,£) plays a very important role as the phase function. The stationary points of 
the phase function 0(k,£) equal to ±y/^ and hence they are real because £ < 0. The signature 
table of the function 

Im 9(k, = (12 Re 2 k - 4 Im 2 k + 12f ) Im k 

depictured on the figure [2j 

Thus Im.9(k, £) > (Im < 0) for k lying in the exterior (interior) of the hyperbola 
3Re 2 fc-Im 2 fc + 3£ = Oof the upper half-plane and in the interior (exterior) of the same hyperbola 
of the lower half-plane. For £ < Im #(&;,£) is negative along (0,ic) and positive along (0, — ic). 
Therefore the jump matrices J(£, t, k) are unbounded (in t) when k G (—ic, ic). Hence we have to 
use the modified nonlinear steepest descent method, suggested in [2], [3], |TT] , and find a new phase 
function g c {k,£), instead of the function 0(k,£), which transforms the original Riemann-Hilbert 
problem to the model RH problem of the finite-gap type. New ^-function leads to the finite-gap 
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Im e(k)=0 



+ 



Re k 



Figure 2: The signature table of the function lm8(k) 

model problems of zero genus for £ < = — c 2 /2 and genus one for = c 2 /3. They are 

explicitly solved by using elementary functions in the first region and the elliptic theta functions 
in the second region, respectively. 

In the region £ < £_ we shall use the following g-function: 

g c (k, = Q(k) + 12£X(k) = (4k 2 - 2c 2 + 12£)X(k) (3.6) 

where X(k) = \/(k — ic) (k + ic) = \Jk 2 + c 2 . This function has the asymptotic behavior similar 
to the phase function 9(k,£), i.e. 

g(k,£) =4k 3 + I2£k + 0(k- 1 ), k ->■ oo. 

The differential of this function can be written in the form: 

. ,, 12k? + (6c 2 + 120*. , (k + X)k(k-X),, 

dg ^ () = — m — * = — m — *■ 

where, evidently, A(£) = \/— £ — c 2 /2. In order to define the boundary value £_ we take into 
account that the phase function g c (k,£) is acceptable until zeroes — A(£) and A(£) are different. 
When they coincide (are equal to zero ) and become complex conjugated then g c (k,£) does not 
work and new phase function must be introduced. Hence we have to put A(£_) = 0. Then 
= — c 2 /2, that gives x < —6c 2 t, and the phase function g c (k,£) will be useful for oo < £ < 
In what follows very important role plays a signature table of the function Img c (/c,£) for 
different values of £. Borderlines between different domains are described by equations: 

&2 = 0; k\ = and |k 2 | < c; 

r 2 \ / r 2 

.2 ;„2 , e , C \ / ;„2 ;„2 - - C 



which are equivalent to lmg c (k,£) = 0. The signature table of the function Img c (k,£) can 
be obtain by using for example "MAPLE" and it is qualitatively depicted on the figure [3] for 
— oo < £ < and on the figure H] for £ = 

Thus g c (k,£) as the function on fc has the following properties: 
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Figure 3: The signature table of the function lmg c (k,£) (£ < 



g c (k,£) is analytic in C\[— ic, ic] 
g c (k,O = 0(k,O + O(k- 1 ), k^oo 
9t(k,£)+g-(k,£) = 0, ke(-ic,ic) 



• 9c(k,£.) = 0(y/k =F ic), k — >■ ±ic. 

The Riemann-Hilbert problem for the matrix M(£, t, k) with the jump contour S = MU [ic, — ic] 
have to be considered now with new (13. 6 j) phase function g c {k,£,)- Let us define new matrix function 

M«(e,t,A;) = M(e,t,fc)G«(e,t,fc), 

where G (1) (£,t, fc) = e^^-^'^K The function MW solve the following R-H problem: 

M w (£,<,*) = M| 1} (£, t , k) J« (£,t,k), M« (£,t,k)^I, k^oo 

where 



' / i _ r (fc) e -2«9e(fe,e) \ 

\ -r(fc)e 2i ' 5c < fc '*> 1 + |r(A;)| 2 J 



jW(£,t,k) = < 



k e 



( 



k e (0,ic) 



\ f(k)e 2it9 ^ 1 / 



/ 1 f(k)e- 2lt ^ k & \ 



I V 



k E (-ic, 0) 




Figure 4: The signature table of the function lmg c (k,£) (£ = £_) 



Further we would like to transfer the jump contour from the real axis. To do so we use the 
following factorizations of the jump matrix on the real axis: 




1 —r(k)e 



-2itg c (k,£) 







1+ \r(k)\ 2 

l + |r(A;)| 2 y 



\ / 1 

—r(k)e 2it9c ^ k '^ 



(3.7) 

\ 

1 



V l + r{k)r(k) I 



(3. 



It is easy to see, the first (second) factor in the first line (13. 7p is decrease as t — > oo in the domains 
where Im g c (k,£) > (Im g c (k,£) < 0). In the second line (13. 8 j) the first(third) factor is decrease 
as t — > oo in the domains where Im g c {k,£) < (Im g c {k,£) > 0). To remove the diagonal terms 
in the second factorization we use a diagonal transformation: 

where some analytic in C\[— A(£), A(£)] function 5(k,£) must be defined. Then the jump matrix 
J {2) (£,t, k) for k G [-A(0, A(0] takes the form 



■/ (2) (£,U) 



'l ± r-r^ 1 (-±(1 



1 + |r| 
1 
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2\-l 



(1 + |r| 2 ) 



\ / 1 °\ 

_ re 2itg c (k£) 



\(l + |r| 2 )<£ ) 



If we chose the function 6(k,£) in the form: 



fc + A(0 



(3.9) 



where 



X(A;,0 = exp 



j_ > Vl + [r(d)l 2 

2™ y 

-A(0 



s — k 



and 



z, = i-ln(l + |r(A(0)| 2 ) 



(3.10) 



(3.11) 



then 5 + = 8-(l + \r\ 2 ) and, therefore, the middle matrix factor become trivial. Thus the jump 
matrix J^(^, t, fc) has a lower/upper factorization for k G [— A(£), A(£)] and a upper/lower factor- 
ization for fc ^ [— A(£), A(£)]: 



J (2) (£,U)= < 



'l a(A;)6(A;)^(jb^)e- 2i ^( fc .O\ / l o\ 

1 J \-a(k)b(k)5l 2 (k,Oe 2itgcik '° 1/ 

fcG[-A(0,A(0] 



1 0W1 r(k)5 2 {k)e- 2it ^ k '0 

-r(k)5- 2 {k,0e 2it9 ^ 1 lo 1 



fc£[-A(0,A(0], 



where we use the identity: 



r(fc) 



l + | r (fc)|2= fl (W*)- 
For fc G [— ic, ic] the jump matrix J^(^,t, k) takes the form: 







f(fy e >t(gl{k,£)+g- c (k,0) e it(g+(k,0-9c (k,0) J ' 
e it(gt(k,€)-gc(k,0) 



k G [0,ic], 
k G [-ic,0\. 



To remove unbounded in t exponentials we have to put 

g;(k,0+9t(k,0 = 0, k G Hc,ic]. (3.12) 
Chosen above function g c (h,£) satisfies this condition and the jump matrix takes the form 

-2itg+(k,0 







J (2 \^k)- [ f{k)5 - 2{k) e 2itgUk,0 I 



k G [ic, 0] 



-2itg+(k,0 f(k)5 2 (k) 



k G [0,ic] 
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Q 7 




Figure 5: The contour S^ 3 ) for the M^(^, t, /c)-problem 

on the contour [ic, — ic]. The jump contour Y,^ for M^(^, t, /c)-problem is the initial one, i.e. 
£( 2 ) = £. 

Let us define a decomposition of the complex fc-plane into eight domains fli, fl 8 separated 
by their common boundary S^ 3 ^ as it is shown on the figure [5j The contours L\, Lq (L3, L4) 
range from the point A(£) (— A(£)) to infinity along the rays argfc = ±7r/4 (argfc = 7r =F 7r/4); 
the contours L2 and L5, passing through the points ±A(£), form a closed oval containing segment 
[—ic, ic] . Then the next transformation is: 



M^{i,t,k) = M^\i,t,k)G^\i,t,k), 



where 



G^\i,t,k) = { 



1 

'1 o N 

1^ 

1 -r(A;)<5 2 (fc,Oe- 2it9c{fc ' 5)N 



6 Oi U fi 3 



5- 



(3.13) 



1 a(fc)6(fc)5 2 (fc,Oe- 2i * 9c(fc ' f)N 



I a(A;)6(A;)r 2 (A;,Oe 2i * 9c(fc ' f) 1 



, k eVL~, 



k e fi 8 . 



(3.14) 
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The G( 2 )-transformation gives to the following RH problem: 

Mi 3) (£,*,*) = Mf (£, t, k) (£,£, k), M { y(Z,t,k)=I + 0(k- 1 ), k -> 00 

on the contour depicted on the Figure |5] with jump matrices J( 3 '(£,t,k) which are equal to 
the identity matrix on the real axis, they coincide with matrices G^ 2 \k) from fl3.13p - fl3.14p written 
for the contours k E Lj (J = 1, 2, 6). It is easy to see that J^(£, t,k) = I + 0(e~ et ) as t — > 00 
and k G Lj with exception of some neighborhoods of the stationary points ±A(£). Therefore the 
main contribution to the asymptotics comes from the jump matrix on the segment [ic, —ic], where 
it takes the form: 

j( 3 ) &t,k) = G<? " (k) fat, k)G m (k) = 



1 -a + (k)b + (k)5 2 (k,0e- 2[t9t{k ^\ , 2)f Jl a^(k)b^{k)5 2 {k,C)e 2it9t ^ 

1 ) ^ k '[ 1 

k e [0, ic] 

1 °)j^(etk)( 1 °' 

-a + (k)b + {k)5- 2 {k,0e 2it9t{k ^ l) \a_{k)b_(k)5- 2 (k,0e- 2it9t{k '° \ 

k e [0, -ic]. 

Using equalities 1 — f(k)a + (k)b + (k) = 0, 1 — f(k)a + (k)b + (k) = and a + (k)b + (k) = —a-{k)b-{k), 
which follow from the definition of the function f(k), we obtain 



J (3) (£,t,fc) 



-f-\k)5 2 (k,0 
J(k)6-\k,0 

f(k)5 2 (k,zy 

(\-f-\k)5- 2 (k,0 



k e [0, ic] 



k e [0, -ic]. 



Further we would like to obtain the RH problem with constant in k jump matrix. To do so let 
us use the factorization 



F+\k,0 \/0 A (F„(k,$ 
F + (k,0) V V V F~\k,0 



which takes place if F_(k,QF + (k,£) = -if{k)5~ 2 {k,i) for k e [0,ic] and F_(k,g)F + (k,£) = 
ir\k,06- 2 {k,0 ^ot k e[0,-ic]. 

Thus we come to the scalar Riemann-Hilbert problem: find a scalar function F(k,£) such that 

• F(k,£) is analytic outside the contour [ic, —ic] which is oriented from ic to —ic; 

• F(k,£) does not vanish; 

• F(k,£) satisfies the jump relation: 

F_(k,0F+(k,0 = h(k)5- 2 (k,0, k e [ic,-ic], 

where 

h(k) = f-i/W = aZ\k)a- + \k) k e [0,ic], 
[if-^k) = a-(k)a+(k), fce[0,-ic]. 



13 



• F(k,£) is bounded at infinity. 
To solve this RH problem let us put 



'Faux (k, £), k G C+\[0, ic] 



a(k) 

a(k)F aux (k,0, keC-\[-ic,0] 



and use the function X{k) = Vk 2 + c 2 . Since 



and 



we have 



Faux{ki 




"log 


Fauxik-i £) 




X(k) _ 


+ 




X{k) _ 




log F avx (k,Z)~ 




log F aux {k^) 


X{k) 




f 


I X{k) 





log5~ 2 {k) 
X + {k) ' 

loga 2 (k) 



k G [ic, —ic 



keR, 



F aU x(k,£) = exp 



X{k) 
2vri 



loga 2 (s) ds 
s-k X{s) 



X{k) ' 

c 

s-k X+(s) 



Let us note, that F(oo,£) = 1. 

Indeed, as o(oo) = 1 then F(k, oo) = F aux {k, oo) 



^a«x(oo,0 = exp { — 



loga 2 ( S )^y- | logr 2 ( S ,0 



ds 



X + (s) 



As 



• for s G M a(s) > & a(-s) = a(s) 

• for s G (-ic,ic) 5(s,£)>0 & 5(-s, £) = l/5(s, f) 

• for s G (—ic,ic) X + (—s) = X + (s) 



then 



IC 

I loga 2 ( S )^y = and J log 5~\s,0 



ds 



X + (s) 



= 



Finally, F(oo,£) = 1. 
Now, since 

where 



J (3 \£,t, k) = F-^(k,0J mod FZ 3 (k,0, k G [ic, 



-2C 



J 



mod 



i 

1 



the next step is as follows: 



M^\i,t,k) = MW(t,t,k)F-°>(k,t), 
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Then we have 
where 



r(4) 



I keR, 

J mod ke{ic,-ic), 

l + 0(e- £t ) keLj, j = 1,2,. ..,6. 



The analysis of the parametrix solutions near the end points ic, —ic and the stationary points 
±A(£) are very similar to the analysis done in [2] and [13], respectively. In the first case, since 
the local representation of g c (k,£) at the points ic and — ic is characterized by a square root type 
behavior: 



9c(k, ~ g {ic, £)Vk-ic, k ->■ 



ic\ 



9c(k, ~ go{~ic, £)Vk + ic, k -)• -ic, 



the relevant model Riemann-Hilbert problems are solvable in terms of the Bessel functions while 
in the second case of the real stationary points ±A(£) they are solvable in terms of the parabolic 
cylinder functions. The contribution to the asymptotics has the order 0(t -3 / 2 ) ?, 0(e~ ev/ *) ? in 
the first case and 0{t~ 1 / 2 ) in the second case. Therefore we have: 



M^(C,t,k) 



J + °(^)) M(m ° d) (fc), 



(3.18) 



where M^ rnod \k) solves the zero-gap model problem (cf. |17j): 



M 



(mod) 



(k) = M ( ™ od) {k)J imod) k e {ic, -ic), M {mod \k) = 1 + O^k- 1 ), k^oo. 



with constant jump matrix: 



J 



mod 



i 

1 



To solve the model problem let us use the function 

x{k) = 



k — ic 
k + ic 



introduced in the first section. Since x^(k) = ix + {k) on the cut {ic, —ic) the explicit solution of 
the model problem takes the form: 



M {mod) {k) 



K{k) + 



1 



(*) 



x{k) 



x{k) - 
x{k) + 



A 



(k) 



X{k) " VV ' K{k)J 

Finally we have the following chain of transformations of the RH problem: 

M{i,t, k) = M {l \i,t, k)e lt ^- 9c(k ^ , 
M^ l \i,t,k) = M<n(Z,t,k)6°*(k,£), 
M^{t,t,k) = M^tMG^&ttk)]- 1 

M^{x,t,k) = M^{U,k)F^{k,0, 
M {i \x,t, k) = M {mod) {k){I + 0{r 1/2 )). 
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Let us emphasize that any matrix M^\^,t,k) (j = 1,2,...) defines the same functions q(x,t) 
since all bordering matrices are diagonal at the point k = oo. By the Theorem 2.1 q(x,t) = 
2i lim [kM(x/12t,t, k)] 12 . If we denote lim [kM®(x/12t, t, k)] 12 = m[ j 2 \x, t) then, take into ac- 

k— >oo ' k— >co 

count the chain of our transformations and using the equalities m^lf d (x,t) = c/2i, F(oc,£) = 1, 
we have: 

q(x,t) = 2imi 2 (x,t) = 2im ( [2 (x,t) = 2im^ 2 (x, t) 
= 2img ) (x,t) + 0(r 1 / 2 ) 
= 2imS 4 2 ) (x,t)F- 2 (oo,0 + 0(r^ 2 ) 
= 2im™ od (x,t)F- 2 (oo,0 + 0(r 1/2 ) 
= c + 0{t- 1 ' 2 ). 

Theorem 3.1. The solution of the IB V problem / II. Ip -z TOj) for t — > oo in the region — oo < x < 
—6c 2 t takes the form: 

q(x,t) = c + 0(r 1/2 ), 



3.2 Elliptic- wave region —6cH < x < Ac 2 t 
3.2.1 The construction of the (^-function 

We need a function g(k,^) with the following properties: 

1. g is analytic against k in C\[— ic, ic] 

2. 3\im(g(k,O-0(k,O)eC 

k — >oo 

3. a set of points where Im g(k,^) = divides the complex plane into four connected open sets 
and contains the set 1R [J[ic, id] [J[—ic, —id], where d e (0,c) is some number. 

We will find such a function in the form: 

k 

fU C\ f ^{S 2 + {1 2 ){S 2 + d 2 )d S r 

J w(s) 

ic 

where the positive on the real axis 1R function w(s) is analytic in C\[— ic, ic] and numbers fi and d 
have to be determined as functions on £. The integration contour is taken to have no intersection 
with the segment (— ic, ic). It is easy to see that g{k,^) G 1R if k e [c, d] + or k G [c, cf]_. To 
satisfy the requirement g{k,^) G M if k G [— c, — d] + or k G [—c,—d]_ we have to choose the 

ic 

numbers /i and d in such a way that J dg(k) = 0. Besides this requirement makes g{k,!;) being 

— ic 

a single- valued function on C \ [—ic, ic] against k. The mentioned requirement can be written as 
follows: 




(3.19) 



o 
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Let us define the function 



F(ji, d) = J (/i 2 - \ 2 d 2 ) 



1-A 2 

c 2 - X 2 d 2 



dX. 



(3.20) 



It is easy to see that there exists a function // = p,(d) such that F(fj,(d), d) = and < //(d) < d. 
Moreover, //(d) is strictly increasing when d G [0, c]. Indeed, one can check that F(/i, d) is strictly 
increasing against // and is strictly decreasing against d when < // < d < c. Now if < d\ < 
c?2 < c then F(//(di), di) = = F(//(cf2), d-z) < F (/j,(d2) , d\) , that is F(//(di),di) < F(fi(d2),di) 
and, hence, //(cii) < ^{d-i). Furthermore fi(d) is continuous function. It follows from 13. 191 



Ad) 



\ 2 d 2 



1-A 2 

c 2 - A 2 d 2 



dX 



1-A 2 

c 2 - A 2 d 2 



dX 



(3.21) 



o ' o 

Now we want to satisfy the requirement (2) 3 lim (g(k,£) — 9(k,£)) G C. For this enough that 



k— >oo 



dg(k, - d6(k, f) = 0(£r 2 ) as k -> cx). Since d0(Jfe, = 12(fc 2 + and d(/(jfe, = [12(fc 2 + 



c 2 cT 

/x 2 — — h — ) + 0(/c 2 )]dk as — >■ oo, we require that 

2 +{ = ^ + T 



(3.22) 



c ^2 

Equation (13.211) yields that /x(0) = and //(c) = — =. Hence // 2 (<i) + — is vary over the segment 

V 3 2 

5c 2 / c 2 c 2 \ 

[0, -— ] when d is vary over the segment [0,c]. So for any £ G I — — , — I there exists a single 

6 \ 2 3 / 

d = G (0, c) such that 



(3.23) 



It is easy to see from f !3.23|) that <i(£) is continuous as the function on £ G ^— — , — ) . Thus the 
function g(k,£) is completely defined. 



c 2 c 2 



1C 



id 



-id 



Figure 6: The signature table for Im g 



The signature table of the imaginary part of the function g(k,£) is depicted on figure [6] 
The function g(k,£) has the following additional properties: 
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2.a lim (g(k,Z)-e(k,Z)) = 

It is follows from the facts that this limit exists and 

g(k,0-6(k,Z)em. , k G (ic, too) 

and 

g(k,S)-8(k,S)eR , keR 

4. g_(k,£) + g + (k,£) = 0, k e (ic,id)\J(-ic,-id); 

5. g-(k,£)-g+(k,£) = B g (£) > 0, ke (-id, id), where 



IC c 

B g {i)=2 J dg{k,i) = 2A J 



id d 



cr — 



24 / 



dt. 



3.2.2 Transition from initial R-H problem to the model one. 

Now we introduce the chain of the transformations of the R-H problem. 
Step 1 

First, we change the phase function 9(k,£) with g(k,£): 



Then 

and 

where 



M {l \i,t, k) = M(£,t, k)e u{9{k '^ e ^ ))a:i . 

M (1) -> I, k ->■ oo 

M { }\i,t,k)=M^\i,t,k)J^\i,t,k), 



J (1 \£,t,k) = [ ], ke 

_ r ( k yitg(h,Z) l + | r (jfe)|2 



e itB g {£) q 

y^ e it(2 5+ (fc,o+B g (0) e~ itB °® 

e itB g (0 j ( k y-it(2g+ (k,0+B g (£)) 
p -itB a (i) 



k E (0, id) 



k E (0, -id), 
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/(*) 



2itg+(k,£) 



k G (ic, id), 



-2itg+{k,es 







/(*) 



o 2itg+(k,0 



k G (—ic, —id). 



Step 2 

The function g(k,£) and its imaginary part (see signature table of the lmg(k, £) on the figure ((5D) 
suggest a choice of new contour S 2 for R-H problem. Let L± = {k : k = s ± — oo < s < 00} 
be the strait lines. Then S 2 = K U [— ic, zc] U L + U L_. We use the lower-upper factorization (13.71) 
of the jump matrix on the real axis and apply the transformation 



f / 



M (2 » (e, t , fc) = M« (e, t, fc)G (2) k), G* a > (£,t, k) 



0\ 



k G fix 



\ - r (A;)e 2i * 9(fc ' ?) 1 / 



/ 1 r(k)e~ 2U9 ^ \ 



fc G 2 



1 / 



J, fcGfi 3 U^4 



(3.24) 



where domains ttj, j = 1, 2, 3, 4 are indicated on the figure (JZj) and get the new R-H problem: 

M^\i,t,k) = M^(i,t,k)J^(i,t,k), (3.25) 
M {2 \^,t,k) ->■ J, fc ->■ 

with the following jump matrices: 

1 N 



J (2) (£,U) 



-r(A;)e 2i * 9(fc) 1 



fc G L+, 
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'1 f(k)e~ 2it9 ^ s 



.0 1 

' e -2it g+ (fc,£) Q 



/(*) e 

-2itg+(k,0 



/(*) 



Q e 2itg+(fc,0 



k e L_, 
fc G (—id, id) 
fc G (ic, id), 
G (—ic, —id). 




Figure 7: The contour E2 for the R-H problem M' 2 '(£, t, fc) 
We have just used the properties 1-5 of the function g(k,^), the jump relation 

/(Jfe) = r _(fc) -r+(fc), fce(0,ic) 

and symmetry relation 

r(k) = -Kl), fc G (C\(-ic, ic)) |J ((-ic, ic) + (J(- 



— zc, zc _ 



Step 3 

First we note that the function /(&) has the following analytic continuation from the interval 

(—ic, ic): 

f(k) = f + (k), ke(-ic,ic), (3.26) 



where 



2 _L r .2 



\k) 



x(k) 



k — ic 
k + ic 



(3.27) 
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Then we factorize the jump matrix on (ic,id) \J(— ic, —id) as follows: 

/ 1 0\ / 1 0\ 



J^ 2 \i,t,k) ==F7°>(k,Z) 



\f+{k)F%{k,t) ) 



0\ ( 1 



F?(k,£), (3.28) 



\f4k)Fl(k,0 ) 



for k G (ic, id) and 



J^(U,k)=F~^(k,0 



\0 



/+(*) 
1 



i 

1 



/ 1 -Fl(k } Oe~ 2itg - {k '°\ 



J 



\0 



/-(*) 
1 



/ 



for G (— ic, —id). Direct calculations show that it is possible if 

• F(k,£) is analytic outside the contour [ic, —ic] which is oriented from ic to — ic. 

• F(k, £) does not vanish; 

• F(k, £) satisfies the jump relation: 



FZ 3 (k,0 
(3.29) 



F_(k,Z)F + (k,£) 



-if{k) = (a-(k)a+(k)) 1 , k G (ic,id) 
a-(k)a+(k), k G (— ic, — id), 



f(k) 



and 



F_(Jfe,0 = F+(k,C)h(k), fee (-id, id), 
where /i(/c) is some function, which has to be determined. 

F(fc, £) is bounded at infinity. 



ic 





Q 3 L< 




^ Q^)l_7 


Q 3 
















Q 2 




0' 






Rek 

Q 2 


Q 4 L< 




+ Q^Ls 





Figure 8: The contour E3 for the R-H problem M^(^,t, k) 
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To solve this RH problem let us put 

1 

I a 



"Faux 

(k, £), k G C+\[0, ic] 



a(k) 

a(k)F aux (k,£), keC-\[-ic,0] 



(3.30) 



. Since 



F aU x-(k,0 F aux+{k,0 = 1, k G (ic, id){J(-ic, -id) 

Faux + (k,0 = F aux _(k,Oa 2 (k), A; G R 
We use the function w(k) = a/ (k 2 + c 2 )(k 2 + d 2 ). Since 

= 0, k E (ic,id)\^J(—ic, —id), 



~\ogF aux (k,t)~ 




~\ogF aux (k,t)~ 


w(k) 


+ 


w(k) 



and 



~logF aux (k,£)~ 




"log F aux {k, £)~ 


log a 2 (/c) 


w(k) 


+ 


w(/c) 


w(fc) 



fc G R, 



we have that one of the functions which satisfy the last two equations is 



2iri J s — k w(s) 



(3.31) 



where 



^(A;,0=e ifeA (l + O(-) ),A;^oo 



A = J_ [ log a 2 (s)ds 
2tt y w(s) 

— oo 



But F aux (k,£) has an essential singularity in infinity. To remove this singularity we introduce an 
Abelian integral of the second kind 



oo 

ft(k) = d + J 



{s 2 + e )ds 
w(s) 



so that 



Q(k) = k + O ( - ) , A;^oo 



id 



I 



(s 2 + eo)ds 
w(s) 



As 



fi(ib) - fc e R, fc g 
fi(fc) - fc g «R, fc G 



(3.32) 
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then do = and 



Then 



(s 2 + eo)ds 



wis 



f2 + + f2_ = on (ic, id) \^J(—ic, — id) 
fi_ — f2 + = B^i < on (id, —id) 



Define 



F aux (k,£) = exp 



w(fc) f In a(s)ds 



<A(on(*,o, 



7ri J (s — k)w(s) 

— oo 

Then F(k,£) defined by (13.301) has the following additional properties: 

• lim F(k,£) = 1, fc ->■ oo 

• F_(ife,0 = F+(fc,Oe _iSn( * )A($) , fc G (-id, id), 

By using the factorizations (I3.29p . (I3.28j) we get the following RH problem: 

(£,t,k) = (£, t, k)G® (£,t,k), Mi 3) (£, t, k) = Mf (£,t,k) J (3) (£, t, fc) , 

oo 



where 



/ i F 2 (fc,Qe- 2 ^^^ \ 



G?®(^t,A;) = F-^(A;,0 



/(*) 
\ 1 
/ 1 °\ 

,2itg(k,£) 



/ 



fc G fi 6 [J f2 s 



and 



/ i F 2 (fc,e)e- 2if g( fe ^ \ 



^ S) (e,*,*) 



/(*) 

\ o 1 

/ 1 o\ 



k E L 



7, 



, k E L$ 



\ f(k)F*{k,t) ) 



\0 1 

/ 1 o\ 

2itg(k,£) , fc G L e 

1 



/ 



V f(k)F*(k,0 ) 
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/ 



—r(k)e 2itglyk '^ 



/ 



/ 



, k e Li 



\ F*(k,0 

e (i« 9 (0+iB n (0A(0)a 3> he (-id, id) 



-r(k)F 2 (k,Q 

e 2itg(k,0 



\ 
/ 



k G L 2 



(3.40) 



( • o)' ke(ic,id)\J(-id,-ic) (3.41) 



Step 4 

Now we introduce a model problem jVfi mo< ^ = jVf| mo< ^ j( mod ) ^ M( mod ) / as k oo, where 

' e (i« 9 (0+iB n (OA(0)<7 3> fc e (_ id> id ) 

(3.42) 



j( mod )(e,t,A;) 



i 

1 



, k G (ic, id) U( — id, ~^ c ) 



To solve it we need a notion of the Riemann surface X, which is induced by 

w 2 (k) = (k 2 + d 2 )(k 2 + d 2 ), 

with cuts along (ic,id) and (—id, —ic). On the first sheet of this surface w(0) > 0. We introduce 
the a-cycle and the b-cycle as it is shown at the picture: 

The basic of the holomorphic differential forms is given by the differential form 



bj = 2ni- 



dk 
w(k) 
f dk 



Then 



dk 
w(k) 



T = T 



(0 == 



w(k) 



2m 



dk 
w(k) 



< 



We introduce the Poincare theta function: 

Q(z) = Q(z,r(0)= e W {K(i)m 2 + zm\ 

It has the following property: 

S(z + 2nin + t(Z)1, r(0) = 0(z, r(£)) exp j"^)/ 2 - *z} , 
Then we introduce the Abel mapping on the surface X: 

p 

A:X^C/ (2niZ + t(£)Z) A(P) = J 



n G N, I G N 



(3.43) 
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Now we introduce the functions ip(., £), £) : {the first sheet of the X} — > C 

B(A(k) - A(Dj) -K- itB g {i) - zS n A(0) 



^,0 



8(-A(fc) - AjDj) -K- itBgjS) - i£ n A(fl) 
e(-A(k) - A(Dj) - K) 



J = 1,2, 



(3.44) 



where .Di = (0, —cd) is on the second sheet and D\ = (0, cd) is on the first sheet. A(D\) = —A(D2). 
K is the Riemann constant of the surface X. The integration contour in ( I3.44p is taken on the 
first sheet and not to intersect the interval (-ic,ic). These functions have the following properties: 



ip_(k,0 = if + (k,0e itB ^ )+lB " A ^ 
V>-(fc,0 = ip+(k^)e~ itB ^~ iBnA ^ 



k G (ic, id) [_J(— id, —ic) 



, k G {id, —id) 



1 1 

k — zc\ 4 f k + id\ 4 



Define a function 

7(fc)=7(fc,0 = 

Then the solution of the model problem of the 4 can be produce as following: 

MiT° d) (U,k) Mir d) (U,k) 



k — id 



k + ic 



M {mod) (£,t, k) 



Vi\oo,0 _ 

^(oojf) 



Let us note that 



7(^0 + 



M (mod) (£,t, jfe) 
1 



7(^,0 



7(^,0- 



7(^,0- 



7(^,0 



A(oo) 



71 i 

Y 



7(^0 



7(^,0 + 



7(^,0 



(3.45) 
(3.46) 

(3.47) 



7T2 r 

if = h - 

2 2 



A{0 



± 



T Til 

T 2 + T 



Taking it in (13.471) we have: 

M {mod \£,t,k) 



( M (mod) ^ ^ k ^ ^ ( M (mod) ^ ^ k ^ 
( M (mod) t; fc j j ^ ( M ( mod ) (£, t, fc)) 



(3.48) 



25 



and 



(t + \) *MP) - y - if5, - iB n A) m 



11 2 „ . . . „. 711 



(M {mod \^k)) 



7 - - 

7 



e(A(p)--) e(rtB 9 + s n A) 

2 



7T? 

Q(-A(P)---itB g -iB n A) Q(0) 



2 Q (-A(p)--) e(fts fl + s n A) 

2 



(M( morf )(e,t,A;)) 21 = 



1 » 771 



1--J Q(A(P) + -- itB g - iB n A) Q(0) 
~ e(A(P) + -) e(itB fl + 5 n A) 

1 \ -7T-7 

7 + e(-A(P) + y- itB g - iB n A) Q(0) 
22 ~ 2 0(-A(P) + y) e(itB g + B a A) 



We can also solve the model problem by using a function 



Then we have that 



M {mod \^t,k) 



(^f (mod) ^ t? k ^ ^ ( M {mod) t? k ^ 



(M^(e,t,fc)) u = 



1 \ 7 - i 

a + aJ e (A(p)- y -itg g -»g n A) eM 

2 



( M (™«0(e,t,fc)) 12 = 



A - - 



77% 



2 0(_A(P) + -) e(itB g + BnA) 

2 

Ay (^( p ) + y " itB 9 ~ e(™) 
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0(A(P) - -) QiltB 9 + 

2 



(M^(e,t,A;)) 22 = 



A + l)e(-A(p) + ^-»tg g -»g n A) eM 

e(_A(p) - -) Q{ltB v + B ^ A ) 

2 
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Comparing fl3^ . fl330|) . fl33T]) .f l332|) with f l3~54j) . fl3~55]) . fl3~56]) . (l3~5Tj) respectively we take that 



2 v y A(fc) ,v y 7(fc) 



©(O) c + d 

eR = V^ (3 ' 59) 

Then 

q mod (x, t) := lim 2zA: (^-, t, k) - i) = lim 2ik t, k) - i) = 

fc->oo \ \ LZt J J 21 fc-s>oo V V iZt J J 12 

8(7ri + ^ g (Q + ^ n (QA(0) e(vrz) 
1 J 6(0) Q(itB g + iB n A) 

9(ttz + ttB g (Q + ^ n (QA(Q, r) 6(tu + iU, r) 

0(rtS fl (O+i5 n (OA(O,r) VC e(iC/,r) 

where U = tB g + B n A. 

Qmod can be expressed in terms of Jacobi elliptic functions: 

q mod (x, t, k) = (c + d)dn [K{m) ( - + 1 ) \m) = . C ~ x (3.61) 



(3.60) 



where 



v v dn ( Kirn) — \m 

V n 

-2irK(l-m) , s 

r = -A— L me 0,1 

2 

K W = I V, T^TTT (3.62) 

1 — msm [v)) 2 



o 



3.2.3 Asimptotics solitons 



c 2 c 2 
As £ — >■ — — then cf — >■ and r — > — oo. Then q m od(x, t) — > c as £ — )■ — — . 

C 2 

If we direct £ to — then r tend to and theta-functions in (13.601) are slowly- convergent. But 
3 



we can use the Poisson summation formula and rewrite the formula f)3.60p in rapidly convergent 
theta-functions. 

e(z,r) = e [——) \j— ( ex P^rJ ( 3 - 63 ) 

Then 

q mod (x, t) = v^^exp ( -f - + ^ ( - + 1 J J V prjj \ J ( 3.0 I ) 
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4tt 2 

where r* = 

T 

Now we have to expand r*, -B 9 , £>q, A in some series which are convergent when £ is tend to 

3 ' 



7 w(fc) y 

6 d 



V(c 2 -y 2 )(l/ 2 -rf 2 ) TT V 2 16 



Let us introduce the method of solving the inverse Jacobi problem (??). Let X be a Riemann 
surface of genus g, {%}j =1 , {6j}| =1 -homological basis of surface X, {ujj}? =1 is the basis of homolog- 
ical differential of the surface X such that J u>i — 2ni5ji and B-is matrix of b-periods of the basis 

a 3 

of holomorphic differentials, A : X — > l-niL 9 + Bl* 9 -is Abelian mapping. Let / be a meromorphic 
function on the Riemann surface X with poles {Qj}™ = \i D = {P/}f =1 -is nonspecial divisor which 
doesn't contain any pole of /. We know that Theta(A(P) — A(D) — K, r) is Abelian integral of 
the third kind with zeroes in points of divisor D. 

Let us integrate a differential form f (P)d In Q(A(P) — A(D) — K,r) along the border of the 
fundamental polygon of the Riemann surface X. 

Then we get 



E/( P i) = 2^E / f(Ppi-i2 res QAf( P )d\n&(A(P) - A(D) - K,r)) 

1=1 1=1 I j=l 



(3.66) 



There are exist curves ju- 



c (x - 4 A) _ V + Q / .o g ( log,A as v ;= l _ 3J ^ Q (3>6?) 



& 8 

such that if point (x, t) lies between two lines •y n (2n-5) and r y 7T (2n+5) then 

q m od{x,t) = O (v 1 - 5 ) as v^O < <5 < 1 
and if point (x, t) lies between two lines 7 7r (2ri+i-<5) and 7 7r (2 n +i-<5) then 

q mod (x, t) = 2 ° (1 + O (v + v^-V)) , < 5 < 1 
cosh(arg) v v 



where 

arg = 2c(x - AcH) + (2n + 1) log(32c 3 t) + ~^ ^ + 27ricA (c) + 

log 

. c(x-4c 2 t) ^ /log 2 (-logMV. , 
+ (2n+l)log 1 > +o( s A ^ — ) , 



*(s) 



log (u) 
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3.3 Vanishing dispersive asymptotics (x > 4c 2 t) 



To study asymptotic behavior of the Riemann-Hilbert problem RH xt in the region x > u 2 t we 
have used well-known technics from [T2], [32], [2]. The large time asymptotics of the solution in 

1/1 k 
4\k + !f 2 

stationary points of the phase function 9(k) are real and equal to ±£. We have 



this region is defined by the phase function 6{k) = 7(7 + 7^)) where £ 2 = t/Ax. Indeed, the 



\k\ 2 - £ 2 

Img(fc)= lmk 

Therefore lm9(k) > (Im8(k) < 0) for k lying in the lower (upper) half-disk and out of the upper 
(lower) half-disk defined by the circle \k\ 2 = £ 2 (Figure ??). For £ 2 < \E\ 2 = l/Au 2 (that is, for 
x > co 2 t) and for k G 7 U 7 the jump matrix J^'(x,t, k) tends to the identity matrix as t — > 00. 
Hence the contour 7 U 7 does not contribute to the main term of the asymptotics, which is defined 
by the stationary points ±£ and has the order 0(t _1//2 ). This asymptotics of the solution was 
done in [?]: 

Theorem 3.2. The solution of the IBV problem (??)-(??) for t — > 00 in the region x > u 2 t has 
a quasi-linear dispersive character, i.e. it is described by the Zakharov - Manakov type formulas: 



q(x, t) =2y - exp ^2iVxi - log Vxi + i<p(£)} + 

+ 2a/ ^ *^~^ exp {-2iVxt + irj(-£) log Vxi + iy?(-g)} + o(r 1/2 ), t -> 00, 
where the functions rj(k) and <p(k) are given by the equations 

v(k) = ^\og(i- P 2 (k)), e 

(f(k) = 3r](k) log 2 — argp(A;) — argT(— ir](k)) H — / log \s — k\dlog[l — p 2 (s)]. 

4 ' 7T J_ £ 



t 



Here T(—ir)(k)) is the Euler gamma-function, and p(k) 



K 2 {k) - 1 



K 2 {k) + 1 

3.4 Zakharov— Manakov region x > 4c 2 £ 

We have the following chain of the transformations of the R-H problem (12.4(1 : First we use 

M« (£, t, k) = M(£, t, k)GM(£, t, k), M W (£, t, k) = M| x) (£, t, k) J (1) (£, t, fc), (3.68) 

1 

GW(£,t,k) = I I , k 6 fli (3.69) 



1 V(^)e- 2ite ^ 



kett 2 (3.70) 



1 
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7, k e n 3 (J fi 4 



(3.71) 



J (1) (U,k)=\ , keL, 

\ _ r ( k yite(k,t) j / 




(3.72) 




1 - r (k)e- 2lte ^ \ 



(3.73) 



= 7, fc e (— ic, ic) 

Now we consider a model problem: Mi mod) (£, t, k) = M [ ™ d \£, t, k)J (mod \£, t, k), 

M^ mod \i,t,k) ->■ 7, A;^oo 



(3.74) 



where 



j( mod \Z,t,k)=I, ke(-ic,ic) 



(3.75) 



M (mod) ig trivia u y solvable: 



M (mod) (£,£,£) = 7 
gmod (x,t) := lim 2tk(M^ 
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Remark 3.1. 7n i/ie paper [1] the asymptotic behavior of the solution of the problem (11), 
(11) was studied in a neighborhood of the leading edge x = u 2 t in terms of asymptotic solitons. 
The problem of matching of the elliptic wave with the asymptotic solitons and these solitons with 
the vanishing self-similar wave is much more complicated and will be considered elsewhere. 
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